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ABSTRACT 


Numerical  machine  calculations  have  been  made  that  determine 
the  effect  of  space  charge,  composed  of  mobile  charge  carriers, 
on  the  drift  and  diffusion  of  those  carriers  in  an  electric  field. 
An  appai'tencltt'.dif fusion  of  carriers  arises  from  the  difference  equa¬ 
tions  used  to  solve  the  transport  problem,  even  though  no  diffusion 
term  is  included  in  the  transport  equations.  This  calculated  dif¬ 
fusion  is  shown  to  be  exact  in  the  absence  of  space  charge  and  in 
the  limit  of  infinitely  small  calculation  intervals.  Numerical  re¬ 
sults  illustrating  the  effect  of  space  charge  upon  apparent  diffu¬ 
sion  coefficients  and  drift  velocities  are  given  for  a  fictitious 
gas.  Calculations  show  good  agreement  with  experimental  data 
taken  in  two  other  laboratories  to  measure  attachment  coefficients 
in  gases.  Exploratory  calculations  have  indicated  the  usefulness 
of  this  model  for  calculations  in  solid-state  transport . 

1 .  INTRODUCTION 


One  method  of  measuring  the  mobility  of  electrons  in  a  gas 
subject  to  an  applied  field  utilizes  short  pulses  of  ultraviolet 
light  incident  upon  the  cathode .  Successive  pulses  of  electrons 
leave  the  cathode  as  a  result  of  photons  created  mainly  near  the 
anode  in  the  discharge  and  the  photoelectric  effect  at  the  cathode. 
The  mobility  then  may  be  deduced  from  the  time  between  successive 
current  peaks.  In  the  course  of  calculations  simulating  these  ex¬ 
perimental  conditions,  it  was  found  that  the  time  between  calculated 
current  peaks  was  considerably  greater  than  the  electron  transit r 
time  obtained  directly  from  the  mobility  used.  In  seeking  the  ex¬ 
planation  of  this  xmexpected  result,  it  was  discovered  that  the  dif¬ 
ference  equations  used  in  the  calculations  introduced  a  "diffusion" 
of  charge  carriers  down  the  steep  concentration  gradient  that  is 
created  in  the  gap.  This  apparent  diffusion  appeared  despite  the 
fact  that  the  diffusion  terms  had  been  omitted  in  the  formulation  of 
the  problem. 

This  report  will  show  that,  in  the  absence  of  ionization  and  of 
space  charge,  the  difference  equations  that  were  used  actually  give 
a  correct  solution  of  the  diffusion  equation  in  the  limit  of  infin¬ 
itely  small  calculation  intervals.  Illustrative  calculations  will 
show  the  effect  of  space  charge  due  to  injection  of  a  large  pulse  of 
carriers  upon  the  apparent  mobility  and  diffusion  of  the  pulse.  Com¬ 
parisons  with  experimental  measurements  used  to  determine  attachment 
coefficients  in  gases  show  the  usefulness  of  the  calculations.  It  is 
suggested  that  the  calculation  technique  will  also  be  of  value  in 
studying  transport  properties  in  the  solid  state. 

2.  MATHEMATICAL  TECHNIQUES 

The  one-dimensional  continuity  equations  governing  the  transport 
of  two  types  of  carriers  in  a  material  are 
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where 


and 


J_  =  p_p  E  +  eD  &n  /3x 

J  =  p  p  E  -  eD  . 

+  +  4-  Sn+/9x 


In  these  equations  J  and  p  =  ne  are  the  current  density  and  the  charge 
density  of  the  carriers,  respectively,  p  is  the  mobility,  and  D  the  dif¬ 
fusion  coefficient  of  the  charge  carriers,  and  E  is  the  electric  field 
intensity.  The  negative  and  positive  sign  subscripts  refer  to  carriers 
of  indicated  sign  of  charge.  In  the  absence  of  space  charge,  the  well 
known  solution  to  these  equations  for  an  injected  pulse  of  NQ  carriers 
per  square  centimeter  at  the  appropriate  electrode  (x  =  0)nat'  tS=!'0  is 


n(x,t) 


N 


s/4TTDt 


exp 


(x-pEt)3 

4Dt 


The  formulation  used  to  compute  the  temporal  growth  of  current  in 
a  gas  has  been  presented  previously  (ref  1) .  Detailed  consideration  has 
also  been  given  to  the  errors,  the  stability,  and  convergence  properties 
of  the  difference  equations  used  for  machine  computation  (ref  2) .  In 
brief,  (1)  and  (2)  with  added  ionization  terms  but  neglecting  diffusion 
are  solved  together  with  Poisson's  equation 


dE 

dx 


n  ) 


and  suitable  initial  and  boundary  conditions. 

Omitting  the  ionization  terms,  the  difference  equations  used  for 
the  numerical  computations  are 


n  (x,  ,t  +  At)  -  n  (x  ,t  ) 
—  k 7  n  —  v- 3  n 


k7  n 


At 


n-(xk7V  _  n-(xk  ~ 


Ax 


and 


n  (x  ,t  +  At)  -  n  (x  ,t  ) 
4-  k7  n _  +  k  n 

At 


n  (x,  +  Ax,t  )  -  n  (x, ,t  ) 
+  k _ 7  n  +  k7  n 

Ax 


Where 

x^  =  kAx 5  k=0,  1,  2,  ...,  M 
t  =  nAt;  n=0,  1,  2,  3,  ... 


Her©  the  number  of  calculation  intervals  across  the  gap  of  distance 
d  is-  I,  i.e.,  d  =  l^x,  and  the  drift  velocities  are  given  by  equa- 
tions  of  the  form  v  =  |jE.  The  difference  equations  were  chosen 
unsymmetrlcal  with  respect  to  x  for  reasons  of  stability  (ref  2). 
The  time  interval  is  chosen  by  the  use  of  the  input  parameter,  F, 
and  the  equation 


At  =  F  - r - -7 — rvT  O) 

max  fv_(x,t),  v+(x,t.)J 

where  stability  conditions  require  F  to  be  less  than  unity  and  v_ 
and  v  to  be  non-negative  (ref  2),  We  will  assume  the  usual  case 
that  v_  >  v  and  set  p  =  F,  q  =  1  -  F,  p*  =  F*  =  |_l  F/u  and  q*  =  1  -  p* . 
Then  in  the+absence  of  space  charge,  (7)  and  (8)  becomes,  using  (9) 

n  (k,  n  +1)  iaf  pn  (k  -1,  n)  +  qn  (k,  n)  (10) 

n  (k,  n  +  1)  =  p*n  (k  +  1,  n)  +  q*n  (k,  n)  (11) 

+  +  "  + 


Each  equation  describes  the  intact  movement  of  the  distribution 
through  a  distance  given  by  the  fraction,  p  or  p*,  of  one  distance 
grid  point  in  the  appropriate  direction  for  each  time  step.  Owing 
to  the  similarity  of  these  two  equations,  only  one  needs  to  be  con¬ 
sidered,  and  (10)  is  chosen. 

Before  the  rigorous  mathematical  derivation,  it  will  be  helpful 
to  go  through  the  actual  steps  the  computer  would  use  for  a  simple 
problem.  The  initial  condition  is  illustrated  in  the  first  diagram 
in  figure  1,  labeled  n  =  0,  corresponding  to  t  =  0.  The  density 
N(k,  n)  is  unity  for  k  =  1  and  zero  elsewhere,  as  shown  by  the  solid 
line.  For  F  (or  p)  =  1/2,  this  triangular  distribution  moves  intact 
one-half  grid  space,  as  shown  by  the  dashed  curve.  However,  the  com¬ 
puter  only  records  data  at  each  grid  point,  hence  the  density  is  1/2 
at  grid  points  k  =1  and  k  =  2  and  zero  elsewhere  :  .  This  trapezoid 
is  then  the  initial  condition  for  the  second  time  step,  shown  in  the 
second  diagram,  labeled  n  =  1.  This  process  is  repeated  for  each 
time  step.  The  fit  to  the  normal  distribution  becomes  better  each 
time  step,  as  shown  by  the  lower  diagrams  where  the  normal  distribu¬ 
tions  are  shown  as  dashed  curves . 

The  mathematical  derivation  may  now  be  continued.  For  initial 


=  0,  anct  n_(k  ^  0,0)  =  0 

lJk  q11  k  =  b  (k;  n,  p)  N'f  (12) 


and  boundary  conditions,  let 

n_(0,0)  =  N^,  n  (-1,  n) 

n-<k'n)  =  No  k T  (n-k)  j 


then 
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where  b  (k;  n,  p)  is  the  binomial  distribution  (ref  3).  Now,  accord¬ 
ing  to  the  De  Moivre-La place  limit  theorem 


lim  W1  x  1  (k  -  np)2 

b(k;  n,p)  =  -  exp - 5— - 

n,  k-f®  \f27rnpq  pq 


(13) 


providing 


lim 


n,  k-»“ 


■  rfe-  —  0,  where  x  =  — - 

nr1*  ’  k  h 


and  h3  =  npq. 

Thd  right-hand  member  of  (13)  is  recognized  as  the  normal  dis¬ 
tribution.  Comparison  of  (12)  with  (5)  shows  that 


npq  =  2Dt/(Ax)2  and  N  =  N/  Ax 


(14) 


Changing  back  to  the  original  notation,  recalling  that 


n  =  t/At,  k  =  x/Ax,  Ak  =  d/M,  p<=  F,  Q  =  1  -  F 


and 


At  =  Fd/Mv 


we  find 

D_  =  (1  -  F)  v'_d/2M 


This  may  be  expressed  also  in  terms  of  the  gap  voltage,  V  =  Ed,  as 


d_  =  (i-g)  v  .  ^  =  (1  ~  v 

pT'  2M  ’  p+  2M 


(15) 


It  is  seen  that  the  diffusion  coefficients  for  the  negative  and 
positive  charge  carriers  cannot  be  adjusted  independently  to  fit  ex¬ 
perimental  values.  For  gases  where  jj.  «  D  /p.  calculated  from 
(15)  is  much  larger  than  D  /p_,  completely  opposite  to  experimental 
observations.  However,  the  computer  program  has  an  option  to  omit 
electron  transport  properties  (ref  2),  (1)  then  being  replaced  by 
the  stationary  ionization  condition.  In  this  case,  the  time  step  is 
determined  solely  by  v  (x,  t)  in  (9).  This  allows  larger  time  steps, 
and  the  positive  ion  diffusion  coefficient  is  determined  by  F  instead 
of  F*,  as  is  the  case  if  p  >  p  . 

The  limit  on  convergence  for  (13)  is  essentially  that- F  cannot 
be  close  to  unity  (the  pulse  being  transported  without  change  of  shape 
if  F  a  1)  or  to  zero  (no  transport).  A  comparison  of  the  maximum  am¬ 
plitude  obtained  by  setting  p  =  q  =  1/2  in  (12)  with  that  for  the  normal 
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distribution  shows  that  the  approximate  error  is  50  n  ^  percent, 
where  n  is  the  number  of  time  calculations. 

It  is  of  course  preferable  to  choose  F  ^  1/2  to  insure  symmetry  of 
diffusion  in  x.  However,  certain  combinations  of  voltages,  mobilities, 
and  diffusion  constants  lead  to  values  of  M  too  small  for  accuracy 
of  calculations  or  too  large  to  keep  the  computer  time  reasonable. 

Values  of  M  from  12  to  100  and  of  F  from  0.01  .to  0.99  have  been  used 
thus  far.  The  skewness  of  the  binomial  distribution  is  given  by  (q-p)/ 
(npq)-1^2,  which  in  the  present  notation  becomes  ( l-2F)/[nF( 1-F) ]^2 ,  show¬ 
ing  that  the  skewness  decreases  inversely  as  the  square  root  of  n,  the 
number  of  calculation  time  steps.  The  possibility  of  simulating  the 
effect  of  back  diffusion  to  an  electrode  by  using  F  =j:  1/2  is  also  sug¬ 
gested. 

3.  ILLUSTRATIVE  EXAMPLES  FOR  A  FICTITIOUS  GAS 


A  number  of  calculations  have  been  made  to  determine  the  effect 
of  space  charge  upon  the  apparent  mobility  and  diffusion  of  a  ficti¬ 
tious  gas  .  This  fictitious  gas  has  p.  p  =  1  x  10s  and  p  p  =  5  x  103 
cm2  torr/v  sec.  Calculations  will  be  reported  for  a  1-cm  gap,  a  gas 
pressure  of  10  torr  and  an  applied  potential  of  100  v.  Values  of  F  and 
M,  and  therefore  of  D  =  (1-F)  p  V/2M,  varied  widely,  D  being  the  dif¬ 
fusion  coefficient  in  the  absence  of  space  charge  as  indicated  by  (15) . 
Figure  2  shows  the  result  of  calculations  for  M  =  20,  F  =  0.5,  n_(0,0)  = 
10°  cm-3  (N0  =  5  x  104  cm-2),  and  no  ions  in  the  gap.  According  to  (15), 
this  calculation  pertains  to  1.25-v  electrons,  and  the  number  of  charge 
carriers  produces  an  initial  field  distortion,  AE/EQ  =  fE(max)  -  E(min)]/ 
E(ave),  of  0.09  percent.  The  fit  to  (5)  using  DQ  from  (15)  is  excellent 
after  only  a  few  time  steps.  The  pronounced  effect  of  space  charge  dis¬ 
tortion  is  shown  in  figure  3  for  the  same  calculation  parameters,  but 
with  n  (0,0)  =  109  cm-3,  the  initial  field  distortion  then  being  90.5 
percent.  An  apparent  diffusion  coefficient,  D*,  defined  from  D*  =  r2/4t, 
where  r  is  half  the  pulse  width  at  1/e  th  of  its  maximum  value,  is  more 
than  eight  times  greater  than  D0 .  The  time  for  the  maximum  collector 
current  (d  =  1)  decreases  roughly  linearly  with  n_(0,0).  For  the  il¬ 
lustrative  example,  the  resulting  apparent  percentage  increase  in  mo¬ 
bility  is  about  1.5  x  10'8  n  (0,0)  percent.  This  is  illustrated  in 
figure  4  for  various  initial  pulse  sizes.  The  resulting  field  distor¬ 
tion  for  the  example  of  figure  3  is  shown  in  figure  5.  The  initial 
densities  that  can  be  used  are  limited  since  the  computer  is  programmed 
to  halt  calculations  if  negative  fields  are  calculated. 


*  The  contrary  increase  of  time  between  current  peaks,  noted  in  the 
introduction  for  calculations  including  ionization,  results  from  the 
steep,  unidirectional  concentration  gradient  of  electrons  present  after 
the  first  gap  traversal. 
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DISTANCE  -CM 

Figure  2.  Transport  of  an  electron  pulse  in  fictitious  gas  across 
a  gap.  Points  indicate  computer  results,  while  solid 
lines  are  plotted  from  equation  (5),  using  equation  (15) 
to  determine  D0.  Space  charge  effect  is  negligible.  The 
times  in  nanoseconds  are  indicated  on  each  curve. 
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Figure  3.  Effect  of  space  charge  upon  transport  in  ficti¬ 
tious  gas.  Conditions  are  identical  with' thoi 
of  figure  2  except  injected  pulse  is  1000  tin* 
larger.  Parameter  is  time  in  nanoseconds 


ANODE  CURRENT  DENSITY,  MICROAMP/CM 


TIME, MICROSEC. 

Figure-  4.  Collection  current  as  a  function  of  time  for 
various  initial  pulse  sizes.  The  parameter 
for  the  various  curves  is  the  value  of  Hq. 


FIELD,  VOL-T/CM 


DISTANCE, CM 

Figure  5.  Field  distortion  across  the  gap  for  the  example  of  figure 
The  parameter  is  again  time  in  nanoseconds. 
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The  fractional  apparent  change,  of  diffusion  coefficient, 

(D*-D  )/D  ,  seems  to  be  a  function  of  the  fractional  field  distor¬ 
tion  ?or  ?hose  calculations  made  so  far.  The  functional  relation¬ 
ship  appears  to  be  approximately 

(D*-D  )/D  =  A(D  .(i)  (AE/E  )  1 ‘SB 

o  o  o  o 

where  A(D  ,  p)  has  varied  from  about  2  to  20  for  each  fraction  ex¬ 
pressed  in  percentage.  Two  examples  illustrating  this  variation 
are  shown  in  figure  6. 

Runs  made  with  n  (0,0)  =  n  (0,0)  at  n  (0,0)  values  used  above 
show  considerably  different  field  distortions  but  the  changes  of 
the  apparent  diffusion  coefficients  are  generally  reduced  only 
about  10  percent  at  most. 

4.  ELECTRON-ATTACHING  GASES 

Under  certain  experimental  conditions,  the  effect  of  diffusion 
and/or  space  charge  may  be  important  in  determining  attachment  co¬ 
efficients,  a,  in  electron-attaching  gases.  Normally,  these  condi¬ 
tions  are  avoided  as  far  as  possible  (ref  4). 

A  number  of  calculations  have  been  made  to  compare  with  unpub¬ 
lished  measurements  of  A,  V.  Phelps  of  the  Westinghouse  Research  Lab¬ 
oratories.  The  experimental  method  has  been  published  (ref  4).  A 
pulse  of  electrons  forms  an  initial  distribution  of  negative  ions 
that  drift  across  a  gap  under  a  constant  applied  field  and  are  col¬ 
lected  to  give  a  temporal  current  pulse.  The  gas  used  was  3  per  cent 
O2  in  He.  The  experimental  conditions  chosen  for  comparison  with  cal¬ 
culations  were  such  that  both  space  charge  and  diffusion  effects  were 
noticeable  and  are  not  typical  of  conditions  used  to  obtain  published 
attachment  coefficients.  They  include  V  =  1.84  v,  p  =  72.4  torr 
(measured  with  gas  at  195°K) ,  d  =  2.54  cm  (i.e.,  E/p  =  0.01  v/cm-torr), 
estimated  area  =  5  cm8,  and  a  temperature  =  195°K,  which  yields  D/p  = 
0.0168  ev.  The  width  of  the  initiating  electron  pulse  was  3  msec  and 
a  collector  duty  cycle  of  6  percent  was  used.  From  data  taken  with 
more  favorable  conditions,  tentative  values  of  v_  =  74  cm/sec  and  <x  = 
0.61  cm-1  were  selected.  A  tentative  initial  distribution  was  chosen 
and  calculations  were  made  with  M  =  28  and  F  =  0.49.  On  the  basis  of 
these  results,  a  better  initial  distribution  was  chosen*  and  new  cal¬ 
culations  made  with  adjusted  a  values.  The  fit  to  the  experimental 
data  with  a  =  0.66  and  v_  =  74  cm/sec,  and  the  initial  distribution 
as  given  in  figure  8,  is  shown  in  figure  7.  The  diffusion  of  the  pulse 
as  it  drifts  across  the  gap  is  also  shown  in  figure  8.  The  initial 
field  distortion  for  the  distribution  used  is  ^  5  percent  as  shown  in 
figure  9,  which  also  shows  the  field  distributions  at  subsequent  times. 


*  Following  the  suggestion  of  Dr.  A.  V.  Phelps,  an  initial  distribution 
was  calculated  by  taking  a  time  average  of  the  computed  density  at  each 
gap  position  over  the  3-msec  duration  of  the  pulse. 
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Figure  6 .  Proportional  change  in  apparent  diffusion  coef f i 
cient  as  a  function  of  the  maximum  field  dis¬ 
tortion. 
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Figure  7. 


Comparison  of  calculated  with  experimental  measurement  of 
collected  negative  ion  pulse  in  oxygen  mixture.  The  ex¬ 
perimental  data  (ref  4)  are  shown  as  solid  curve  and  cal¬ 
culations  as  individual  points .  The  initial  peak  in  the 
experimental  curve  is  due  to  electrons  that  cross  the  gap 
without  attachment. 


DENSITY  OF  IONS 


Figure  8;  Initial  distribution  of  negative  ions  used  for  calcu¬ 
lation  of  figure  7,  together  with  calculated  dis¬ 
tributions  for  later  tines,  indicated  on  each  curve 
in  milliseconds . 
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In  order  to  separate  the  effect  of  diffusion  from  that  of 
space  charge  upon  the  collected  pulse,  the  initial  distribution 
was  reduced . by  a  factor  of  100  for  an  initial  distribution  cor¬ 
responding  to  ct  -  0.58.  Doshring  (ref  5)  has  calculated  the  effect 
of  diffusion  upon  the  collected  pulse  and  gave  an  approximation 
that  the  apparent  attachment  coefficient,  Q /,  is  related  to  the  true 
coefficient  by 

a'  =  a(l  +  QD/pE)  (16) 


For  the  illustrative  case  one  obtains  aJ /CL  =1.013  from  (16), 
using  v  =85  cm/sec.  For  the  space-charge-free  run  on  the  compu¬ 
ter,  a  value  of  a' /a  =  1.017  was  obtained  initially,  a  value  that 
dropped  slowly  to  1.013  at  about  15  msec,  in  substantial  agreement 
with  (16).  For  the  corresponding  run  with  the  densities  near  that 
used  experimentally,  cc'/a  initially  was  equal  to  1.04,  but  dropped 
to  unity  at  about  10  msec.  This  shows  that  at  a  current  density  of 
about  10~ia  amp/cmg,  space  charge  causes  a  greater  error  in  the 
attachment  coefficient  than  does  diffusion  at  195°K. 

The  calculations  bear  out  the  statement  (ref  4)  that  the  true 
transit  time  is  intermediate  between  the  time  of  the  ion  current 
peak  and  the  time  at  which  the  ion  current  has  dropped  to  half  its 
peak  value. 

A  second  experimental  method  was  used  by  JJorrison  and  Edelson 
(ref  6)  to  measure  attachment  coefficients  in.SFg  at  10  torr.  They 
measured  the  total  current  due  to  the  movement  of  ions  within  a  gap 
of  0.8  cm.  Using  intense  illumination  they  obtained  considerable 
distortion  due  to  space  charge.  They  also  presented  a  calculation 
of  the  effect  of  space  charge  distortion,  but  omitted  the  diffusion 
of  carriers.  The  conditions  assumed  in  their  calculations  were  used 
with  the  present  computer  program.  A  comparison  of  the  present  re¬ 
sults  with  their  calculation  and  with  their  experimental  data  is 
shown  in  figure  10.  The  agreement  between  the  two  calculations  is 
satisfactory,  but  shows  the  effect  of  diffusion  on  the  apparent 
mobility.  Better  agreement  with  the  experimental  data  could  be  ob¬ 
tained  with  different  values  of  the  attachment  coefficient  and  mobil¬ 
ity  of  the  ions.  The  calculations  were  made  with  M  =  32,  F  =  0.8 
and  V  =  150  v.  It  is  noted  that  these  values  yield  D  /p  =  0.47  ev, 
which  is  probably  too  high  a  value.  A  better  value  of  the  ion  mobil¬ 
ity  would  undoubtedly  require  a  smaller  diffusion  energy  to  fit  the 
experimental  data.  These  calculations  were  made  before  (15)  was  de¬ 
rived. 

5 .  SOLID-STATE  TRANSPORT 


The  injection  of  minority  carriers  into  a  semiconductor  and 
the  measurement  of  the  pulse  at  a  collector  after  drift  and  diffusion 
is  utilized  to  measure  mobilities  in  a  semiconductor.  Normally 
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pulses  of  several  sizes  are  used  and  the  drift  time  is  extrapolated 
to  zero  pulse  current  to  account  for  space  charge  effect  (ref  7). 

Although  the  present  formulation  v/as  developed  for  gases,  it 
may  be  applied  with  reservation  to  the  solid  state.  The  present 
program  is  limited  to  a  dielectric  constant  of  unity;  hence  carrier 
densities  must  be  multiplied  by  the  dielectric  constant  of  the  solid 
to  relate  to  the  space  charge  effect.  Recombination  and  trapping 
are  omitted  from  the  formulation;  the  former  could  easily  be  intro¬ 
duced  while  the  latter  is  more  complicated.  Further,  in  most  semi¬ 
conductors,  doping  creates  a  large  concentration  of  majority  car¬ 
riers,  which  are' neutralized  by  fixed  charges.  Such  fixed  charges 
could  easily' be  introduced  in  a  computer  program. 

Exploratory  calculations  have  been  made  to  indicate  the  use¬ 
fulness  of  the  present  program  in  determining  the  effect  of  space 
charge  and  diffusion  upon  mobility  measurements  in  germanium.  Calcu¬ 
lations  made  assuming  carriers  of  only  one  charge  show  that  under 
typical  experimental  conditions,  a  space  charge  sufficient  to  broaden 
the  pulse  to  the  observed  width  distorts  the  pulses  to  nearly  rectan¬ 
gular  shapes .  The  use  of  sufficient  carriers  of  the  opposite  sign  to 
neutralize  the  gap  as  a  whole,  but  distributed  uniformly  across  the 
gap,  leads  to  quite  different  field  distortions,  and  to  more  normally 
distributed  pulses.  Calculations  were  made  with:  V  =  40  v,  d  =  2.0  cm, 
p_  =  3800  cm2/v  sec,  p  =  1800  cm2/v  sec,  M  =  100,  F  =  0.87,  and  the 
total  number  of  initial  carriers  of  one  charge  in  the  gap  per  square 
centimeter  varied  up  to  8  x  10s .  The  collected  pulses  for  several  in¬ 
itial  pulse  sizes  are  shown  in  figure  11.  Comparison  with  experimental 
measurements  shows  some  disagreement.  Experimental  measurements  show 
a  pronounced  increase  in  the  time  taken  for  the  pulse  maximum  to  be 
collected  as  the  pulse  size  increases.  The  calculations  made  with 
carriers  of  only  one  charge  show  a  slight  decrease  in  this  time,  al¬ 
though  the  total  time  for  pulse  collection  increases  by  the  correct 
order  of  magnitude.  Calculations  made  with  initial  overall  charge 
neutrality  show  peaks  delayed  in  time.  By  overcompensation  of  the 
injected  carriers,  the  pulse  is  further  delayed  as  shown  in  figure  11. 
This  suggests  that  a  program  utilizing  fixed  charges  would  enable 
attainment  of  good  agreement  with  experimental  data.  It  is  further 
suggested  that  temporary  trapping  could  account  for  further  delay  be¬ 
cause  carriers  may  be  removed  from  the  early  portion  of  the  collected 
pulse  and  added  to  the  later  portions. 

It  is  believed  that  the  present  model  for  calculations  of  diffu¬ 
sion  effects  in  the  presence  of  space  charge  can  easily  be  adapted 
to  make  valid  calculations  for  semiconductors.  By  proper  selection 
of  a  spatial  distribution  of  fixed  charges,  a  semiconductor  junction 
may  also  be  simulated. 

5.  DISCUSSION 

It  has  been  shown  that  the  present  program  is  very  useful  in 
interpreting  the  effect  of  space  charge  upon  the  transport  properties 
of  charge  carriers. 


TIME  -/-*« SEC 

Figure  11.  Collected  pulses  of  holes  calculated  as  a  function 
of  size  of  injected  pulse  of  holes  in  Ge.  Pulse 
size  is  indicated  by  relative,  not  absolute,  num¬ 
bers.  Calculations  with  only  holes  present  are  shown 
dashed,  while  initially  overall  neutral  calculations 
are  shown  as  a  solid  curve,  and  those  for  an  over¬ 
compensated  electron  distribution  by  the  dot-dashed 
curves.  The  hole  transit  time,  tp,  is  indicated. 
Results  for  diffusion  at  two  temperatures  are  shown 
in  two  cases. 
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It  is  natural  to  ask  why  the  exact  second-order  equations  (3) 
and  (4)  including  diffusion  directly  are  not  used.  The  main  diffi¬ 
culty  is  that  an  additional  boundary  condition  is  needed  for  each 
of  the  two  equations.  Considerable  effort  has  been  made  to  deter¬ 
mine  these  two  extra  boundary  conditions  for  steady-state,  rather 
than  time-dependent ,  calculations.  No  solution  to  the  problem  has 
been  found  that  has  been  considered  worthwhile  to  test  on  the  com¬ 
puter.  It  is  possible  that  the  dynamic  problem-  is  easier,  but  any 
set  of  equations  must  be  investigated  to  insure  that  the  solution 
is  exact  in  the  limit  of  small  differences.  The  present  program 
has  proved  to  be  exact,  including  diffusion,  in  that  limit. 

The  chief  deficiency  of  the  present  program  is  that  the  posi¬ 
tive  ion  temperature  must  be  higher  than  the  electron  (or  negative 
ion)  temperature'  as  shown  by  (15) .  Just  the  reverse  is  the  usual 
case  in  gases.  In  many  cases,  large  concentration  gradients  apply 
to  the  carriers  of  one  sign  only,  in  which  case  the  diffusion  of 
carriers  of  the  other  sign  is  not  important  anyway.  However,  this 
fundamental  deficiency  of  the  present  calculations  indicates  that 
the  second-order  equations  should  be  solved  numerically. 

It  is  planned  to  publish  at  a  later  date  the  results  of  calcu¬ 
lations  where  the  ionization  coefficients  are  not  set  equal  to  zero. 
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